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Chapter  1 

Second-Order  Nonlinear  Optics  and  its  Applications  in  Electro-Optic  Modulators 

1.1  Theoretical  Background 

1.1.1  Photonics,  Electronics  and  the  motivation  for  EO  polymer  research 

For  several  years  the  world  has  relied  on  electronics  for  information  processing. 
However,  we  are  nearing  the  limits  of  traditional  electronic  signal  processing  and 
transmission.  Optical  technology  will  be  required  in  order  to  gain  processing  rates  above 
50  GHz  and  a  transmission  distance  greater  than  1  meter.  The  optical  technology  is 
described  as  photonics,  which  is  the  analog  of  electronics.  This  means  that  photons, 
rather  than  electrons,  will  be  used  to  acquire,  store,  transmit  and  process  information. 
Many  photonic  applications  are  already  in  use.  These  include  fiber-optic  communication 
lines,  diode  laser  sources,  organic  liquid-crystal  displays  and  optical  memories.1 

There  are  many  advantages  of  photonics  over  electronics.  Optical  switching,  for 
example,  occurs  on  a  femtosecond  time  scale.  No  electronic  process  can  come  close  to 
matching  that  time,  but  optical  processes  can.  And  with  optical  signals,  there  is  no 
electrical  and  magnetic  interference,  so  there  is  less  crosstalk  between  channels. 

In  the  past  few  years  there  has  been  a  demand  for  more  bandwidth.  The 
explosion  of  multimedia  applications  through  the  Internet  brought  a  requirement  for 
higher  bandwidth  in  order  to  run  fast  and  smooth.  Glass  fibers  are  most  commonly  used 
to  transmit  data  at  high  speeds  because  they  have  an  inherent  large  bandwidth  of  more 
than  lOTHz.  However,  a  decrease  in  the  rate  of  transfer  occurs  at  the  electric-to-optic 
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signal  conversion  intersection  between  fibers.  The  data  gets  bottlenecked  in  these  areas.2 
Electro-optic  modulators  are  a  combination  of  electronics  and  photonics  that  help 
alleviate  the  bottleneck. 

These  modulators  serve  important  fiber-optic  and  wireless  functions,  including: 
satellite  receiver  systems  collecting  high-quality  rf  signals  that  have  been  distributed 
long  distances,  connecting  remote  cellular  radio  systems  with  high-quality  rf  signals  over 
long  distances,  and  local  area  networks  where  low-distortion  rf  signals  are  needed. 
Electro-optic  modulators  also  find  uses  in  other  areas,  such  as  frequency  mixing,  light 
modulation,  optical  logic,  millimeter  wave  signal  generation,  optical  beam  steering,  rf 
detection,  phase  control,  power  splitting,  and  wavelength  division  multiplexing 
(WDM)4’5 

1.1.2  Inorganic  Crystals  vs  Organic  Polymers 

The  field  of  nonlinear  optics  (NLO)  took  off  with  the  coming  of  the  laser.  Laser 
beams  have  strong  oscillating  electrical  fields  that  are  capable  of  polarizing  a  medium 
nonlinearly  corresponding  to  the  field  strength.  The  medium  can  then  serve  as  a  new 
optical  field  source  with  different  frequencies  and  modified  propagation  qualities. 

All  forms  of  matter  could  exhibit  nonlinear  optical  effects  if  the  incoming  optical 
field  were  strong  enough.  However,  with  today’s  laser  sources,  only  a  few  materials 
display  significant  nonlinear  optical  effects. 

There  are  basically  two  types  of  nonlinear  optical  materials,  inorganic  crystals 
and  organic  polymers,  and  classification  depends  on  the  origin  of  nonlinearity.  The 
inorganic-crystal  species,  or  ferroelectric  crystals,  include  GaAs,  KTiOPCL  (KTP), 
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LiNb03  and  LiTa03.  The  nonlinear  property  of  these  materials  comes  from  electrons 
shared  by  the  bulk  material,  such  as  in  metals  and  semiconductors.  The  electronic 
characteristics  of  the  bulk  medium  determine  the  optical  nonlinearity.  The  polymer 
materials  include  organic  crystals  and  polymers.  The  nonlinearity  in  these  materials 
comes  from  the  molecular  structure  and  geometrical  arrangement  of  the  molecules  in  the 
condensed  medium. 

Nalwa  and  Miyata  list  ten  characteristics  ideal  nonlinear  optical  materials  should 

have6: 

1 .  Large  nonlinear  figure  of  merit  for  frequency  conversion 

2.  High  laser  damage  threshold 

3.  Fast  optical  response  time 

4.  Wide  phase  matchable  angle 

5.  Architectural  flexibility  for  molecular  design  and  morphology 

6.  Ability  to  process  into  crystals,  thin  films 

7.  Optical  transparency  (no  absorption  at  fundamental  and  SH  wavelengths) 

8.  Ease  of  fabrication 

9.  Nontoxicity  and  good  environmental  stability 

10.  High  mechanical  strength  and  thermal  stability 

Polymer  materials  have  a  few  advantages  over  inorganic  crystalline  materials. 
Organic  molecules  can  be  synthesized,  processed  and  integrated  in  ways  inorganic 
crystals  cannot  be.  This  synthetic  freedom  allows  researchers  to  optimize  nonlinear 
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responses.  Molecular  materials  also  exhibit  exceptional  mechanical  properties,  thermal 
and  chemical  stability,  optical  transparency  and  high  laser  damage  thresholds. 

Inorganic  crystals  are  based  on  positive  and  negative  ions.  The  ionic  lattice 
deforms  relatively  slowly  in  response  to  laser  light  interaction.  This  deformation  leads  to 
heat  dissipation,  which  limits  the  cycle  time  of  these  devices.  These  resonant  optical 
nonlinearities  have  a  response  time  established  by  the  lifetime  of  the  excitation.  On  the 
other  hand,  polymer  materials  have  large,  nonresonant  optical  susceptibilities  derived 
from  ultrafast  excitations  and  a  large  displacement  of  conjugated  ^-electrons.  The 
nonlinearity  is  limited  only  by  the  width  of  the  driving  laser  pulse. 

Another  advantage  is  that  polymer  materials  have  low  dc  dielectric  constants 
around  the  order  of  3  as  compared  to  28  for  LiNb03,  for  example.  This  results  in  a 
minimization  of  phase  mismatch  between  optical  and  electrical  pulses  allowing  a  large 
operating  bandwidth  modulation.5  Electro-optic  modulators  have  nearly  the  same  phase 
velocities  at  optical  and  microwave  frequencies.  A  few  years  ago,  waveguide 
modulators  were  proven  to  operate  at  frequencies  as  high  as  150  GHz. 

Other  advantages  include  the  ability  of  polymer  materials  to  be  deposited  directly 
onto  a  number  of  substrates,  even  conforming  to  flexible  substrates,  which  is  not 
practical  with  LiNb03.  Combine  this  with  the  ability  to  operate  modulators  at  very  low 
voltages  (<1  V)  compared  to  LiNb03  modulators  (~5V).  All  this  adds  up  to  why  so  much 
research  and  effort  has  been  put  forth  in  the  search  for  high-performance  polymeric 
electro-optic  materials. 

1.2  Introduction  and  Theory  of  Nonlinear  Optics 
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1.2.1  Microscopic  Nonlinearity 

A  linear  response  is  simply  the  lowest-order  approximation  describing  a  process. 
A  process  is  termed  nonlinear  if  the  response  to  an  input  physically  changes  the  process 
itself.  Nonlinear  optics  deals  with  the  interactions  between  light  waves  and  matter  or 
light  waves  and  other  light  waves.  When  electromagnetic  radiation  intermingles  with  a 
nonlinear  optical  material,  the  optical  characteristics  of  the  material  change,  which 
transforms  the  properties  of  the  interacting  light  itself. 

To  describe  nonlinear  optical  phenomena  requires  looking  at  the  origin  of 
electronic  polarization.  At  the  atomic  level,  the  electromagnetic  field  E  (E  =  coscot) 
exerts  a  force  F  on  an  electron: 

F  =  -eE  (1.1) 

The  original  state  is  perturbed  from  its  equilibrium  position,  and  the  electron  density  is 
displaced.  This  displacement  gives  rise  to  a  charge  separation  or  electron  density 
polarization  P: 

P  =  -er  (1.2) 

where  r  is  the  displacement  of  charge.  The  polarization  is  driven  by  the  applied  electric 
field,  so  it  is  easiest  to  picture  the  displaced  system  as  a  Lorenz  Oscillator.7  The  electron 
acts  as  a  harmonic  oscillator  bound  to  a  nucleus  by  a  spring  with  natural  frequency  coo- 
Next,  consider  only  the  linear  response  in  which  the  restoring  force  from  the  nucleus  is 
approximated  as  linearly  proportional  to  the  displacement.  The  equation  of  motion  is 
then  written  as  a  second  order  differential  equation: 
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„  2  d2r  dr 

eE  -  ma>o  r  =  m  — —  +  21  — 
dt 2  dt 


(1.3) 


The  terms  on  the  left  side  of  the  equation  represent  the  Coulomb  force  and  the  restoring 
force.  The  right  side  contains  a  damping  term  with  coefficient  T  to  account  for  energy 
dissipation  during  the  polarization.  Solving  the  equation  for  the  displacement  r  gives: 

r  =  -— — - - - rE„e"  +  c.c.  (1.4) 

m  a>Q  -  2iT  o)-o) 

Equation  (1.4)  reveals  sinusoidal  behavior  in  time  with  increasing  displacement 
as  the  frequency  of  the  field  approaches  the  natural  frequency  of  the  oscillator.  The 
atomic  polarization  is  then  written  as: 


P=  - - - -  E0eicot  +  c.c.  (1.5) 

m  co0  -2  iToo-co 

This  equation  shows  two  things.  First,  the  molecular  polarization  occurs 
sinusoidally  with  the  identical  frequency  ©  to  the  driven  field.  Second,  the  amplitude  of 
polarization  is  proportional  to  the  strength  of  the  applied  field.  Rearranging  equation 
(1.5)  into  its  real  and  imaginary  parts  gives: 


?Re  : 


2  2 
0)q  -cd 


m(o)l-co)  +4T  co 


E0eio>t  +  c.c. 


(1.6) 


Plm  “ 


2  Tco 


Eoe'®*  +  c.c. 


m  {col  -  cor)  +4 To 


(1.7) 


The  real  part,  or  in-phase  part,  of  the  polarizability  shows  that  a  new  source  of  radiation 
with  the  same  frequency  as  the  fundamental  wave  evolves  as  the  electromagnetic  wave 
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propagates  through  the  medium.  The  imaginary  part  reflects  the  out-of-phase 
polarizability  and  is  attributed  to  the  damping  effect.  This  happens  when  a  charge 
doesn’t  respond  immediately  to  the  applied  field  creating  a  phase  difference  between  the 
polarization  and  the  electric  field  results.  Refractive  index  and  birefringence  are  linear 
optical  effects  that  result  from  linear  polarization.  Electron  polarization,  which  occurs  on 
the  femtosecond  time  scale,  is  the  fastest  of  all  polarizations.  For  this  reason  organic 
nonlinearity  is  superior  to  inorganic  crystals  in  terms  of  speed. 

By  similar  reasoning  the  equation  of  motion  describing  nonlinear  optical 
processes  is  derived  after  adding  an  anharmonic  term  ar2  to  equation  (1.3).  This  term 
acts  as  the  first-degree  correction  to  the  nonlinear  dependence  of  the  restoring  force  on 
the  electron  displacement.  Nearly  all  second  order  NLO  effects  are  ascribed  to  this 
added  term. 

eE  -  ma>o2r  =  m^-j-  +  2T  —  +  ar2  (1.8) 

dt2  dt 

By  adding  the  anharmonic  term,  the  exact  solution  to  equation  (1.8)  can  no  longer 
be  derived.  The  way  around  this  is  to  assume  that  the  anharmonic  contribution  to  the 
total  polarization  is  inconsequential.  Then  the  solution  to  equation  (1.8)  can  be 
approximated  as  a  power  series  in  charge  displacement  r.  Accounting  for  the  first  two 
terms  of  the  expansion: 

r  =  ri+r2  (1.9) 

The  solution  to  rj  is  practically  identical  to  that  of  equation  (1.3).  Solving  for  the 
second  term  requires  approximating  ar2  as  arj.  Then  equation  (1.8)  becomes  a  linear 


differential  equation,  which  leads  to  an  exact  solution  in  terms  of  a  component  at 
frequency  2co  and  another  at  zero  frequency  or  dc: 
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r2(0)  =  -2a 


Eja )  1 

cOq  (g)q  -  2iT co- co2) 


(1.12) 


The  r2(2co)  term  shows  in  the  exponential  numerator  that  a  new  source  of  sinusoidal 
radiation  at  twice  the  original  frequency  co  is  being  generated  simply  by  including  the 
second  order  correction  term  ar2  for  the  restoring  force.  This  phenomenon  is  called 
second-harmonic  generation  (SHG).  It  is  important  to  see  that  the  amplitude  of  the  SHG 
is  proportional  to  the  square  of  the  applied  electric  field  strength.  The  r2(0)  term  shows 
that  along  with  SHG,  a  frequency  independent  polarization  is  concurrently  produced. 

This  phenomenon  is  called  optical  rectification. 

Up  to  this  point  both  linear  and  nonlinear  optical  processes  have  been  defined  as 
the  result  of  electron  density  redistribution  after  interaction  with  light  waves.  The  linear 
effects  are  based  on  an  approximation  that  the  restoring  force  is  linearly  dependent  on  the 
electron  displacement.  In  this  case  the  electron  redistribution  is  symmetric.  When  a 
higher  order  approximation  to  the  displacement  is  made,  as  in  the  case  of  nonlinear 
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optical  processes,  the  electron  density  redistribution  is  asymmetric.  To  get  a  non-zero  p 
the  molecule  must  have  acentric  symmetry. 

A  second  way  to  approach  optical  nonlinear  polarization  is  to  express  it  as  a 
power  series  of  the  field  strength  E: 

p  =  aE  +  pEE  +  yEEE  +  •••  (1.13) 

where  a  is  the  first-order  linear  polarizability,  P  is  the  first  order  molecular 
hyperpolarizability  and  y  is  the  second  order  hyperpolarizability.  Substituting  the 
electric  field  of  a  plane  light  wave,  E  =  E0cos(cot),  into  equation  (1.13)  gives: 


1  9 

p  =  p0  +  aEocos(utf  -kr)+  —  PEq  [1  +  cos(2<s t  -  kr )] 


+  yEcf 


'3  1 

—cos (a>t  -kr)  +  —  cos(3 cot  -  3 kr) 

4  4 


(1.15) 


The  third  and  fourth  terms  reveal  new  radiations  are  generated  at  frequencies  2©  and  3®, 
respectively.  The  third  term  also  shows  a  frequency  independent  dc  polarization. 
Incidentally,  the  3©  frequency  is  called  third  harmonic  generation.  Again,  if  p  is  to  be 
non-zero,  the  molecule  must  be  asymmetric. 

Another  way  to  present  the  NLO  process  is  to  examine  how  two  different  waves 
with  electric  fields  Ei  and  E2  interact  with  NLO  material  electrons.  Using  two  laser 
beams  with  frequencies  ©1  and  ©2,  the  P  in  equation  (1.14)  becomes: 

PEiCos(©it)E2COs(©2t)  (L16) 

Equation  (1.14)  then  looks  like: 

PEiCOS(©it)E2COS(©2t)  =  —  PEiE2C0s[(©i  +  ©2)t]  +  —  pEiE2COs[(©i  -  ©2)t]  (1.17) 

2  2 
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From  this  it  is  easy  to  see  that  polarizations  occur  at  frequencies  (<»i  +  ©2)  and  (©1  -  ©2). 
Radiation  will  be  re-emitted  at  these  frequencies  with  a  dependence  on  the  coefficient  p. 
Notice  that  when  ©1  =  ©2  the  result  is  both  SHG  and  optical  rectification. 

From  the  results  of  the  previous  equations  it  is  seen  that  the  interaction  of  light 
with  NLO  materials  changes  the  subsequent  properties  of  the  light.  Whether  the  process 
is  looked  at  as  electron  distribution  or  wave-wave  interaction,  the  responses  both  display 
optical  nonlinearity.5 

1.2.2  Quantitative  View  of  First  Order  Hyperpolarizability 

When  the  oscillating  electric  field  of  a  light  wave  interacts  with  the  electrons  in 
an  NLO  medium,  it  perturbs  the  electrons  from  their  ground  state  and  causes  oscillating 
currents  in  the  molecule.  The  new  electronic  state  is  now  a  linear  combination  of  the 
ground  state  and  excited  quantum  states  of  the  un-perturbed  system.  The  electric  field  is 
time  dependent,  which  means  that  both  the  number  of  excited  quantum  states  and  the 
electron  density  distribution  will  fluctuate  with  time.  Simply  put,  the  oscillating  electric 
field  generates  a  time-dependent  polarization.  The  fluctuating  electron  density  can  then 
be  described  by  time-dependent  perturbation  theory.8 

The  time-dependent  Schrodinger  equation  is: 

HT=i  h&¥/dt  (1.18) 

The  perturbing  Hamiltonian  is: 

H'  =  -erE(osin©t  (1.19) 

The  overall  Hamiltonian  sums  the  unperturbed  molecular  Hamiltonian  Ho  and  the 


interaction  perturbation  H' : 
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H  =  Ho  +  H'  =  Ho  -  erEtaSinwt  (1-20) 

The  r  term  is  the  dipole  moment  operator  summed  over  all  electronic  coordinates: 

r=X>.  (1.21) 

a 

The  total  molecular  polarization,  dipole  plus  induced  dipole,  is  then  written  as: 

p  =  -e(vF|r|vF}  (1.22) 

The  next  steps  involve  solving  the  Schrodinger  equation,  substituting  the  solution  into 
equation  (1.22)  and  collecting  second  order  terms  for  the  polarization  equation.  The  first 
hyperpolarizability  p  is  then  derived: 

The  p'  term  is  the  dipole  moment  difference  operator  defined  as: 

(«'  M  n)  =  (»'  |p|  n)  -  ( g  |jx|  g)  5„.„  ( 1 .24) 

The  5n'n  term  is  the  Kronecker  delta,  and  p  is  the  dipole  moment  operator.  The  dipole 
moment,  p,  describes  the  extent  of  charge  separation.  The  transition  dipole  moment, 
in'  |p|  nj ,  determines  the  amount  of  mixing  between  the  ground  and  excited  states.  When 

n'  =  n,  5n'n  -  1,  and  the  result  gives  the  difference  between  the  excited  state  dipole 
moment  and  the  ground  state  dipole  moment.  If  n'  *n,  Sn'n  =  0  giving  the  transition 


dipole  moment  from  n  to  n'. 


12 


In  1977  Oudar  used  a  simple  two-state  model  to  approximate  p  from  the 
complete  sum  frequency  generation  perturbation  calculation: 


3e2ft20y(Ap) 

P  2m[lF2  -  (2fao)2  JjF2  -(^co)2] 

where 


(1.25) 


Ap  =  pe-pg  (1 .26) 

In  equation  (1 .25)  p  is  a  function  of  the  laser  fundamental  photon  energy  co,  the 
electron  charge  e  and  mass  m,  the  energy  gap  W,  the  oscillator  strength  f  of  the  charge 
transfer  transition  in  the  molecule,  and  the  change  in  dipole  moment  Ap  resulting  from 
that  transition.  To  maximize  P  large  Ap  and / values  are  required.  Alternatively,  a  small 
energy  gap  W  near  2co  would  also  be  desirable.9 

Marder  et  al  later  reduced  equation  (1.25)  to: 

p  _  (l^ee  /q  27) 

(Mge)2 


where  (pee  -  Pgg)  is  the  difference  between  the  excited-state  and  the  ground  state  dipole 
moments.  The  coefficients  pge  and  AEge  are  the  transition  dipole  moment  and  the  optical 
(HOMO-LUMO)  bandgap,  respectively.10 

1.2.3  Importance  of  Refractive  Index 

The  index  of  refraction  is  the  ratio  of  the  speed  of  light  in  a  vacuum  to  the  speed 
of  light  in  a  specific  medium.  The  refractive  index  of  a  material  is  an  important 
component  for  determining  the  second-order  nonlinear  susceptibility  x(2),  which,  along 


with  the  first  hyperpolarizability  p,  is  needed  to  determine  the  second  order  nonlinear 
optical  properties  of  a  particular  organic  material. 

Analogous  to  the  atomic  polarization,  equation  (1.13),  the  bulk  polarization  is 
written  as: 

p  =  x(1>e  +  x(2)EE  +  x(3)EEE  (1 .28) 

where  %(1),  y(2\  and  %(3)  are  the  macroscopic  susceptibilities.  The  coefficient  yX2)  is  a 
third-rank  tensor  describing  the  second-order  nonlinear  optical  susceptibility  and  is  borne 
only  of  noncentrosymmetric  mediums. 

Optical  radiation  causes  a  disturbance  in  the  electronic  distribution  of  a  molecule, 
which  generates  an  electric  field  opposing  the  applied  field.  The  total  internal  electric 
field  D  converts  to: 

D  =  E  +  4tcP  (1.29) 

where  47iP  is  the  internal  electric  field  generated  by  the  polarization.  The  dielectric 
constant  s  is  defined  as  the  ratio  of  D  to  the  applied  field  E,  or: 


E  +  4nP 

8  =  - 

E 


(1.30) 


Rewrite  equation  (1.27)  as: 

P  =  X<$E  (1.31) 

The  effective  susceptibility,  yeff,  is  a  function  of  the  applied  field  representing  the  total 
bulk  polarization.  Substituting  equation  (1.30)  into  (1.29)  exposes  the  dependence  of  the 
dielectric  constant  s  on  the  polarization: 


e((o)  =  1  +  471  Xe/w) 


(1.32) 
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When  working  at  optical  frequencies  and  in  the  absence  of  dispersion,  the  dielectric 
coefficient  of  an  isotopic  material  equals  the  square  of  the  refractive  index: 

r|2(co)  =  s(co)  =  1  +  4nXeAa)  0  -33) 

If  the  optical  field,  E(co),  is  modulated  by  a  dc  field,  E(0),  the  overall  field  E 
acting  on  the  material  is  the  sum  of  the  two  fields: 

E  =  E(0)  +  E(co)  =  E(0)  +  E0cos(cot  -  kz)  (1 .34) 

Combining  equation  (1.34)  with  (1.28)  gives  (ignoring  higher  order  terms): 

P  =  x(1)[E(0)  +  E0cos(cot  -  kz)]  +  x(2)[E(0)  +  E0cos(oc»t  -  kz)]2 

+  X(3)[E(0)  +  E0cos(oot  -  kz)]3  (1.35) 

Expanding  equation  (1.35)  gives: 

P  =  x(1)E0cos(cot  -  kz)  +  E2(0)  +  2x(2)E(0)E0cos(cDt  -  kz)  +  E02cos(cot  -  kz) 

+  x(3)E3(0)  +  3x(3)E2(0)EoCOs(rot  -  kz)  +  3x(3)E(0)Eo2cos2(cot  -  kz) 

+  x(3)E3ocos3(cot  -  kz)  (1 .36) 

Next,  reduce  cos2  and  cos3  terms.  Then  collect  terms  containing  the  frequency  to: 

P(co)  =  x(1)Eocos(cot  -  kz)  +  2x(2)E(0)Eocos(tot  -  kz)  +  3x(3)E2(0)EoCOs(cot  -  kz) 

+  —  X(3)E3ocos(cot  -  kz)  (1-37) 

4 

Factor  and  rewrite  equation  (1.37)  as: 

P(to)  =  Xe^E(0)cos(cot  -  kz)  (1 .38) 

Substitute  the  factored  terms  of  %e//E  into  equation  (1.33)  to  get: 

ti2(co)  =  1  +  47t[X(1)  +  2x(2)E(0)  +  3x(3)E2(0)  +  ^x(3)Eq2] 


(1.39) 
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Designate  the  first-order  linear  refractive  index  as  fto2,  which  conies  from  the  linear 
polarization  r)o  =  1  +  4n%(l>  to  get: 

r|2(oo)  -  r|o2  =  87ix(2)E(0)  +  127iX<3)E2(0)  +  37t%(3)Eo2  (1-40) 

The  left  side  of  equation  (1.40)  can  be  rewritten  as: 

ft2(©)  -  r|o2  =  (n  +  ftoXft  -  r|0)  »  2rio(ft  -  ft0)  0-41) 

substitution  of  (1.41)  into  (1.40)  gives: 


4-tiy  ^ 

ft(c ))  =  T|0  + - E(0)  + 

fto 


6itx(3) 

fto 


E2(0)  + 


3tix(3) 

2ft  o 


E02 


(1.42) 


It  is  possible  to  show  a  direct  relation  between  the  nonlinear  refractive  index  and  the 
light  intensity.  The  definition  of  light  intensity  in  cgs  units  is: 


Eo2=— I(<o) 
Cft 


(1.43) 


Substitute  (1.42)  into  (1.41)  getting: 


„(.)  -  n.  +  E(0)  +  ^(0)  + 

fto  fto  Cfto 


(1.44) 


Rewrite  equation  (1.43)  as: 

ft(co)  =  fto(co)  +  r)]E(0)  +  r)2(0)E2(0)  +  ft2(co)I(©)  (1-45) 

where  fti  is  the  linear  electro-optic  (LEO)  effect  (the  Pockels  effect),  ft2(0)  is  the 
quadratic  electro-optic  effect  and  ft2(co)  shows  the  direct  relationship  between  the 
nonlinear  refractive  index  and  the  light  intensity  (Kerr  effect).  The  ni  term  given  by: 

4ixX(2) 
fti  =  — ^ — 
fto 


(1.46) 
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is  proportional  to  the  second-order  nonlinear  susceptibility  x(2)-  This  means  that  by 
applying  a  dc  electric  field  to  the  medium,  polarization  will  occur.  This  in  turn  alters  the 
refractive  indices,  which  changes  the  properties  of  the  interacting  light,  such  as  a  phase 


shift  or  a  change  of  direction.11 

The  equation  of  the  phase  shift  of  the  material  is: 

7tr|3rVL 


(1.47) 


and 

V„=4V  (1.48) 

r|  rL 

where  A<j>  is  the  electric  field  dependent  phase  shift  of  light  passing  through  an  electro¬ 
optic  material,  V  is  the  applied  electric  field,  h  is  the  electrode  spacing,  r\  is  the  field- 
independent  refractive  index,  r  is  the  electro-optic  coefficient  of  the  material,  L  is  the 
length  of  the  material,  X  is  the  wavelength  of  light  and  V*  is  the  voltage  required  to 
produce  a  n  phase  shift.  The  phase  shift  at  given  V  and  L  is  linearly  proportional  to  q3r, 
which  is  also  known  as  the  figure  of  merit  (FOM)  for  second-order  NLO  materials  for 
light  modulators.4 


1.2.4  Extensive  View  of  Electro-optic  Effect 

The  linear  electro-optic  effect,  or  Pockels  effect,  occurs  when  a  static  field 
induces  a  field-dependent  change  in  the  index  of  refraction  at  frequency  co.  It  is  directly 
observed  only  in  noncentrosymmetric  crystals.12  It  can  be  thought  of  as  a  three-wave 
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mixing  process,  where  two  waves  at  frequencies  ©1  and  C02  interact  to  produce  a  third 
wave  at  frequency  ©3.  The  second-order  susceptibility  has  the  form  x(2)(-©3,©i,©2).  For 
the  Pockels  effect,  ©2  =  0,  giving  x(2)(-®3,©i,0). 

There  are  many  examples  where  the  LEO  effect  plays  an  important  role  in 
technology,  among  them  being  signal  transduction  and  optical  switching.4,13,14’15 
1.2.5  Exclusive  View  of  Electro-optic  Effect 

Consider  a  light  wave  of  arbitrary  polarization  traveling  through  a  birefringent 
medium.  Under  Maxwell’s  equations,  the  optical  field  is  seen  as  consisting  of  two 
mutually  orthogonal  linearly  polarized  electric  field  components.  Each  component  of  the 
light  wave  is  commonly  subjected  to  different  refractive  indices,  which  leads  to 
birefringence.  The  indicatrix  of  refractive  index,  a  potential  energy  surface  that 
characterizes  the  refractive  index  of  the  traveling  light,  is  used  to  explain  this  occurrence. 

To  derive  the  indicatrix,  begin  with  the  expression  of  the  dielectric  constant  e  in 
an  isotropic  medium.  From  equations  (1.29)  and  (1.30),  s  relates  is  a  second-rank  tensor 
relating  the  electric  field  E  and  the  dielectric  displacement  D.  The  vectors  of  D,  E  and  s 
can  be  expressed  in  matrix  form: 

sn  e12  s13 

S21  S22  S23 

S31  E32  e33 

If  the  medium  is  non-dispersive,  the  dielectric  constant  s  tensor  is  symmetric.  In  the 
principal  axes  system  all  the  off-diagonal  elements  are  zero: 


)(EA 

Ey  (1.49) 

AE  z7 
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8n  0  0  E, 


0  0  833  J  l  e2 


(1.50) 


The  diagonal  elements  are  related  to  the  refractive  indices  by: 


T|x=  Vs 


(1.51) 


by  *\/s  22 


(1.52) 


(1.53) 


Zemike  and  Midwinter  define  the  total  electric  energy  density,  W,  in  the  medium 


W=  — 


1  f  D!  D2V  Di 


—  — —  H - —  -4 - - 

i  tlx  fiy  fiz 


(1.54) 


Simplify  equation  (1.53)  with  the  following  relationship: 


x  =  .  x  . ,  similarly  for  y  and  z 


(1.55) 


to  get: 


x  y  z 

T  +  - T  +  —  ~^ 

fix  fiy  fiz 


(1.56) 


This  is  the  equation  for  an  ellipse.  The  elliptical  axes  are  parallel  to  the  principal 
directions  in  the  medium.  The  incident  light  of  wave-vector  S  has  a  fundamental 
frequency  ca.  The  two  orthogonal  polarization  components  are  parallel  to  (IcosO  + 
KsinG)  and  J,  where  I,  J  and  K  are  unit  vectors  in  the  molecular  x,  y  and  z  directions. 


The  first  orthogonal  component  is  defined  as  r|0,  which  is  independent  of  the  propagation 
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direction,  normal  to  S  and  parallel  to  J.  The  second  component  is  termed  t]c.  It  is 
normal  to  both  J  and  S  and  is  dependent  on  the  angle  of  inclination,  0,  of  the  incident 
beam  with  the  optical  axis.  A  positive  uniaxial  medium  results,  when,  in  equation  (1 .55), 
r|x =  r)y  ^  r|z-  The  refractive-index  ellipsoid  for  this  type  of  medium  is  shown  in  Fig.  1.1. 


z 


Fig.  1.1  Refractive  index  of  uniaxial  medium, 
where  S  is  the  direction  of  light  propagation, 
r|o  is  the  refractive  index  normal  to  S  and 
parallel  to  J,  r|e  is  the  refractive  index  normal 
to  both  S  and  J. 


Applying  an  electric  field  to  the  nonlinear  material  alters  the  polarization 
susceptibility  and  the  refractive  index.  On  the  molecular  scale,  the  electron  density  is 
rearranged  anisotropically.  The  overall  effect  of  the  rearrangement  is  a  change  in  the 
shape  or  rotations  to  the  ellipse.  Either  effect  gives  new  values  for  the  refractive  indices. 
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This  means  that  a  beam  of  polarized  light  will  have  its  plane  of  polarization  rotated 
according  to  the  strength  and  direction  of  the  applied  field  subsequently  changing  the 
velocity  and  possibly  the  direction  of  the  light.  This  is  shown  in  Fig.  1.2. 


a 


b 


c 


Fig.  1 .2  Refractive  index  of  an  ellipsoid  in  (a)  the  absence  of  an  applied  electric 
field,  (b)  a  medium  field  and  (c)  a  strong  field 


Due  to  the  applied  electric  field,  additional  math  terms  are  needed  for  a  complete 
description  of  the  ellipsoid.  Equation  (1.55)  becomes: 


f  1  ' 


^2; 


x2  + 


f  1 ' 


VT) 


y2  + 


<  l  > 


Vh2; 


z2  +  2 


f  i  ' 


Wj 


yz+2 


ii2; 


XZ  + 


f  1  ^ 


^2; 


xy  =  1 


(1.57) 


and 
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r21  r22  r23 

r31  r32  *33 

A(n2f 

r41  r42  r43 

A(n2F 

r51  r52  r53 

4\2V 

r61  r62  r63 

E, 

e2 

e3 


(1.58) 


Let  x,  y  and  z  be  parallel  to  the  principal  dielectric  axes  of  the  medium.  In  the  absence 
of  an  electric  field,  equation  (1.58)  turns  into: 

4(x2+y!)+4z2=1  (1'59> 

n; 


The  expression  relating  the  change  in  the  refractive  index  as  a  function  of  the  applied 
electric  field  to  the  new  and  old  values  of  the  refractive  indices  is: 


for  i  =  1  to  6 


(1.59) 


The  ry  quantity  in  the  contracted  form  of  notation  can  be  replaced  by  rm.  The 

relationship  between  the  contracted  and  uncontracted  suffixes  is: 

ij  ->m 
11  ->  1 
22  — »  2 

33  ->  3  (1.60) 
23  or  32  ->  4 
31  or  13  ->  5 
12  or  21  ->  6 

The  rjj  matrix  is  under  the  same  symmetry  constraints  as  the  second-order  NLO 
coefficient  and  is  directly  related  to  it. 

As  an  example  of  how  an  applied  electric  field  changes  the  refractive  index  of  a 
uniaxial  medium,  consider  dipolar  chromophores  in  an  amorphous  polymer  film.  Under 


the  influence  of  a  strong  electric  field,  the  chromophores  will  align  parallel  to  the  poling 
direction,  which  is  perpendicular  to  the  film.  The  point  group  of  the  poled  film  is  oomm, 
or  Ccov-  The  symmetry  operations  have  an  infinite-fold  rotation  about  the  poling  field 
direction  and  an  infinite  number  of  mirror  planes  including  the  poling  field  direction. 
The  electro-optic  coefficient  tensor  nj  becomes: 


^  Ji  H 


=  V  r-E  = 
Z_i  ij  j 


0 

0 

0 

r42 

0 

0 


E, 

E, 


(1.61) 


There  will  be  two  pairs  of  equivalent  indices  due  to  the  symmetry  of  the  medium.  In  the 
absence  of  absorption,  under  the  Kleinman  symmetry  constraint,  rn  =  r23  =  r42  =  r5i 
giving  two  independent  components  of  the  electro-optic  coefficient  r||  =  r33  and  rj_  =  ro, 
which  are  parallel  and  perpendicular  to  the  polar  axis,  respectively. 

For  a  dc  field  applied  in  the  z  direction: 


1  ] 

fO  . 

f  n 

=  A 

—  = A 

2 

U2J 

! 

In  )2 

ln«  J 

—  ri3Ez 


(1.62) 


and 


44] 

fl] 

2 

-  A 

W) 

3 

r33Ez 


(1.63) 


Now  equation  (1.63)  becomes: 
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In  general  rn  is  smaller  than  r33,  and  in  highly  anisotropic  materials,  ri3  approaches  zero. 

An  alternative  derivation  of  the  light  output  phase  difference  comes  from  the 
reflection  technique16: 

Output  =  AIsin2^  (1-65) 

where  AI  is  the  intensity  contrast,  Imax  -  Imm,  and  't'sp  is  the  phase  difference  between  s- 
(TE)  and  p-  (TM)  polarized  waves.  Note  that  the  incident  beam  is  polarized  45  degrees, 
thus  intensities  for  both  TE  and  TM  are  identical.  The  equation  for  the  phase  difference 
is  written: 

^sP=K+^sP  (\.66) 

where  T'sp0  is  the  phase  retardation  due  to  birefringence  and  84%  is  the  induced  phase 
retardation  due  to  the  Pockels  effect  (or  LEO  effect,  in  this  regard,  Kerr’s  optical  effect  is 
assumed  to  contribute  to  the  modulation  neglectably  in  contrast  to  the  Pockels  effect.) 

9¥sp  -  d¥p  -  d¥s  (1.67) 

The  term  8¥p  is  the  induced  phase  retardation  of  TM  polarization  and  $¥s  is  the 
induced  phase  shift  of  TE  polarization.  The  incident  wave  beam  is  defined  as: 

'P  =  ytiL  (1.68) 

where  t|  is  the  refractive  index,  L  is  the  optical  path  length  and  X  is  the  free  space 
wavelength  of  the  modulated  beam.  Introduce  a  small  phase  difference  to  equation 
(1.68): 
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8T  -b(—r\L 

(1.69) 

^  J 

8W  =  ^-(ri8L  +  L8r\) 

(1.70) 

OA 

When  L  = - ,  where  d  is  the  electrode  distance,  the  derivative  of  L  is: 

cos  a 


5L=2dsma6a  (1.71) 

cos  a 


By  Snell’s  law,  sin  0  =  qg  sin  0g  =  r|  sin  a ,  where  0  is  the  angle  of  the  laser  beam 

irradiating  on  the  glass  surface  from  air  (r|=l),  %  is  the  refractive  index  of  the  glass 
substrate,  0g  is  the  refracting  angle  in  the  glass,  r\  is  the  refractive  index  of  the  polymer 
and  a  is  the  refracting  angle  in  the  polymer  media.  For  fixed  angle  measurements, 

O  =  8sin0  =  8r|sina  =  8(r|sina)  (1-72) 

0  =  r|(Sa)cosa  +  (8r))sina  (1-73) 


~  sina  _ 

8a  = - 8r| 

rjcosa 

Substitute  (1.74)  into  (1.71)  to  get: 

2d  sin  a  C 


(1.74) 


8L  = 


cos2  a 


sma  - 

- 8r| 

^  rj  cos  a  ) 

Combine  equations  (1.75)  and  (1.70)  to  get: 

f  2dsin2  a 


2d  sin 2  a 
rj  cos3  a 


Sr| 


(1.75) 


271 

8^  =  — 
X 


L8r|-r|x 


V 


3  Sri 
rj  cos  a  j 


^  2n 

2d  sin  a  ^ 

=  - 

L - r - 

J  * 

^  cos  a  J 

8r| 


(1.76) 
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Recall  equation  (1 .66),  5^  =  -  8\PS . 


6^  =  — 
sp  X 


2d  sin 2  a 

-L  Z 


V 


cos  a 


|(Srip -SiqJ 


(1.77) 


Substitute  (1.77)  into  (1.65).  The  final  output  intensity  is  then  given  by: 


L  output 


Alsin2 

’5%  n 

2  A, 

^  cos  a  ) 

(1.78) 


1.2.6  Relating  Molecular  Hyperpolarizability  to  Macroscopic  Electro-optic  Activity 

The  equation  for  the  macroscopic  polarization  is  given  by  equation  (1.28): 

p  =  x(De  +  x(2)EE  +  x(3)EEE  (1 .28) 

The  second-order  nonlinear  susceptibility,  %(2\  is  zero  for  centrosymmetric  materials. 

The  expression  relating  macroscopic  y(2)  to  P  is: 

X(2)  =  Np/(co,)/(co2)/(co3)r(Q)  (1.79) 


where  N  is  the  molecular  number  density  and  /( ©;)  are  the  local  field  factors  due  to 
intermolecular  interactions.  The  /( coj)  terms  are  derived  from  the  Onsager  correction 
equation: 


/(©O  - 


fel+2)k) 
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(1.80) 


For  an  isotropic  distribution  of  chromophores,  equals  zero.  An  orientation 
factor,  T(Q),  provides  the  projection  of  p  onto  the  macroscopic  framework  established 
by  the  direction  of  the  applied  poling  field. 
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To  observe  second-order  nonlinear  optical  effects,  externally  applied  dc  fields  are 
required  to  induce  Coov  cylindrical  polar  symmetry.  Under  the  Kleinman  restrictions 
only  two  independent  x(2)  tensor  components  are  nonzero: 

X»  =^P/(®i)/((°2)/(0)3)(cos3e}  (L81) 

and 

XS  =  W(®l)/(®2)/(®3)(COSe)(COs2  0)  (L82) 

The  dc  field  removes  the  orientational  averaging  replacing  T(Q)  with  <  >.  The  0 
term  is  the  angle  between  the  dipole  moment  of  the  nonlinear  molecule  (parallel  to  (3) 
and  the  dc  poling  field,  E.  The  macroscopic  electro-optic  coefficient,  r33,  is  related  to  P 
(expressed  in  esu  units)  by: 

2JVp/(©i  )/(ro2)/(©3)(cos3e) 

r33= - — — -  (1.83) 

(».) 

where  (cos30)  is  the  acentric  ordering  parameter.  The  following  chapter  begins  with 
equation  (1.83)  and  deals  exclusively  with  the  details  of  interactions  between 
chromophores. 

1.3  Electro-optic  Polymeric  Device  Requirements 

A  few  ideas  from  nonlinear  optical  material  theory  are  needed  in  order  to 
fabricate  viable  devices.  The  materials  should  have  high  r33  values,  low  optical  loss, 
good  thermal  and  chemical  stability  during  operating  conditions,  and  the  ability  to 
maintain  polar  alignment. 
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For  practical  reasons,  electro-optic  polymeric  devices  need  to  have  operational 
voltages,  V*,  of  less  than  three  volts.  Equation  (1.83)  shows  that  high-P  chromophores 
are  needed  to  produce  high  r 33  values.  Many  disciplinary  fields,  such  as  optical, 
electrical  and  molecular  engineering,  chemical  theorists  and  synthetic  chemists  combine 
to  create,  test  and  integrate  novel  NLO  materials  into  useful  devices.  Over  the  years 
many  high-P  chromophores  have  been  produced,  but  they  do  not  always  lead  to  high 
macroscopic  r33  values.  One  fundamental  reason  was  provided  by  Dalton  et.  al.,  who 
showed  that  electrostatic  interactions  between  high-P  chromophores  could  sometimes 
limit  the  correlation  between  micro-  and  macroscopic  optical  nonlinearity.4  The 
electrostatic  interaction  produces  chromophore  aggregates  which  do  not  contribute  to  the 
electro-optic  effect  when  the  chromophore  loading  density  reaches  a  certain  level. 
Therefore,  the  molecular  architecture  of  the  chromophores  is  very  important  in  order  to 
produce  high-P  chromophores  that  won’t  clump  together  under  optional  loading 
conditions. 

It  is  important  for  devices  to  have  as  little  optical  loss  as  possible.  Optical  loss  is 
classified  either  as  loss  due  to  absorption  of  light  or  due  to  scattering  of  light.  Loss  by 
absorption  has  two  causes.  The  first  arises  in  the  absorption  tail  of  the  interband 
(HOMO-LUMO)  absorption  and  the  second  is  due  to  absorption  caused  C-H  group 
vibrations.  However,  if  needed,  the  hydrogens  of  the  C-H  groups  can  be  replaced  with 
fluorine,  and  a  large  difference  between  the  operating  wavelength  and  interband 
absorption  maximum  can  be  maintained.  This  means  that  the  proper  chromophore  and 
polymer  design  can  keep  absorption  loss  to  a  minimum.  The  light  should  have  a  total 
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attenuation  of  less  than  1  dB/cm.  This  means  the  chromophores  should  not  absorb  light 
at  operational  wavelengths.  Since  many  high-P  chromophores  have  absorption  maxima 
above  650  nm, 18,19  they  are  unfit  for  device  application. 

Optical  loss  due  to  scattering  is  a  different  problem.  A  major  contributor  is  the 
heterogeneity  of  the  index  of  refraction  in  materials.  This  heterogeneity  can  arise  from 
chromophore/polymer  incompatibility  or  process-induced  losses,  such  as  spin  casting, 
electric  field  poling,  waveguide  fabrication,  and/or  cladding  deposition  of  thin  films. 
Another  area  of  optical  loss  arises  from  impedance  mismatch  between  silica  fibers  and 
polymeric  waveguides.4 

The  thermal  and  chemical  stability  of  chromophores  and  polymers  is  essential  for 
both  the  fabrication  and  operation  of  devices.  There  seems  to  be  a  trade-off  between 
high-P  chromophores  and  high  thermal  stability.  For  instance,  increasing  the  71- 
conjugation  length  gives  higher  P  values  but  lowers  the  decomposition  temperature.  It 
has  been  shown  that  the  thermal  stability  of  electro-optic  activity  can  be  correlated  with 
the  difference  between  the  operating  temperature  and  the  glass  transition  temperature 
(Tg)  of  the  chromophore/polymer  composite.19,20 

Stabilizing  the  polar  alignment  of  the  polymeric  electro-optic  materials  critical 
for  device  operation.  Electro-optic  activity  must  be  stable  for  both  long-term  operation 
at  temperatures  up  to  125  °C  and  short-term  operation  at  temperatures  nearing  200  °C. 
There  are  currently  two  methods  to  achieve  this.  The  first  is  to  utilize  high-Tg  polymers 
in  chromophore/polymer  composite  materials.  The  second  method  involves  covalently 
coupling  the  chromophore  and  polymer,  then  inducing  some  type  of  lattice  hardening 
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near  the  end  of  the  poling  process.  The  chromophores  are  either  tethered  to  the  polymer 
backbone,20’21’22  or  incorporated  directly  into  the  polymer  backbone.23,24  The  purpose  of 
this  study  was  to  develop  high- (3  chromophores  with  improved  thermal  and  chemical 
stability  and  potentially  less  electrostatic  interaction  by  way  of  a  new  donor. 
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Chapter  2 

The  Importance  of  Chromophore-Chromophore  Electrostatic  Interactions  and 
Their  Impact  on  Nonlinear  Optical  Properties 

2.1.  Debye,  London  and  the  Electrostatic  Interaction  Models 

The  work  of  Debye  in  the  late  1930s1  resulted  in  theories  still  used  today  to 
describe  the  behavior  of  noninteracting  dipoles  under  the  influence  of  an  external  electric 
field.  Chromophores  in  a  polymer  matrix  are  viewed  as  independent  particles.  These 
independent  particles  form  the  gas  phase  model  of  macroscopic  optical  nonlinearity, 
which,  until  recently,  was  universally  accepted  as  the  calculation  standard  for  the  pp/Mw 
figure  of  merit.  Here  p  is  the  chromophore  ground-state  dipole  moment,  P  is  the  first 
hyperpolarizability  and  Mw  is  the  molecular  weight  of  the  nonlinear  optically-active 
moiety  (chromophore).  Recently  developed  chromophores,  such  as  FTC  and  CZC  , 
with  high  pp  values  on  the  order  of  10'44  esu  and  r33  values  ranging  from  55  to  84  pm  V' 
*,  have  surpassed  lithium  niobate  (r33  «  30  pm  V'1).  Under  the  gas  phase  model,  the 
macroscopic  electro-optic  activity  was  predicted  to  scale  as  pp/Mw.  However,  with  this 
scaling,  the  bulk  nonlinear  optical  coefficients  from  high-P  chromophores  are  well  below 
the  values  expected  from  the  chromophore  hyperpolarizabilities.  The  gas  phase  model 
appears  to  be  inadequate  for  high-P  chromophores,  which  has  led  to  the  examination  of 
other  theories  exploiting  electrostatic  interactions. 

Chromophores  consist  of  aromatic  and  7t-conjugated  subunits.  As  the 
polarization  of  the  electron  cloud  increases,  the  strength  of  electrostatic  interactions  also 
generally  increases.  Chromophores,  then,  will  have  greater  electrostatic  interactions  than 


saturated  molecules.  The  interactions  between  neighboring  chromophores  alter  the 
electron  density  distribution  in  each  molecule,  which  changes  the  molecule’s  chemical, 
physical  and  photophysical  properties.  Thus,  electrostatic  interactions  must  be  taken  into 
account. 

The  work  of  Fritz  London  in  the  1930s4  led  to  the  ability  to  calculate,  using 
statistical  mechanics,  electrostatic  interaction  contributions.  He  showed  that  the 
attraction  between  molecules  varied  with  distance,  r,  as  r'6.  Molecules  with  small  dipole 
moments  and  polarizabilities  have  fairly  small  electrostatic  interactions.  This  means  that 
other  interactions  are  competing,  such  as  solute-solvent,  solute-solute,  and  H-bonding 
interactions.  For  these  molecules  the  London  forces  are  short-ranged.  But  as  the  dipole 
moments  and  polarizabilities  increase,  the  solute-solute  interactions  dominate,  and  the 
overall  electrostatic  interactions  are  longer  ranged.  It  has  been  shown  recently  that  large 
dipole  moments  and  polarizabilities  of  chromophores  result  in  strong  London  forces.5,6,7 
These  forces  are  largely  responsible  for  the  disparity  between  high-p  chromophores  and 
low  electro-optic  coefficients.  This  chapter  will  show  how  London  theory  is  used  both  to 
define  electro-optic  activity  and  reveal  the  importance  of  chromophore  structure  to 
optical  nonlinearity. 

2.2  Theories  of  Electro-optic  Activity  and  Electrostatic  Interactions 

Begin  with  equation  (1.83),  where  r|e  is  the  extraordinary  refractive  index.  The 

ry  A 

r33  term  is  in  esu  units,  where  1cm  /esu  =  4.19xl0‘  m/V. 

2AT5/(o)I)/(o>2)/(coJ){cosJe) 
r33  = - — 4 -  (l-83) 
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The  acentric  ordering  parameter,  <cos30),  defines  the  competition  between  the 
chromophore/electric  field  interactions  and  the  thermal  randomization  energy.  It  is  also 
possible  to  estimate  the  electro-optic  coefficient  as: 


T33  =  3 


(  V 

1  J\o_ 

VeJ 


(rn) 


(2.1) 


but  this  does  not  always  hold  true  with  electric  field  poling8  and  hardly  ever  holds  true 
for  laser-assisted  poling. 

To  determine  (cos30),  start  with  the  Gibbs  distribution  function: 

exp[-t/(Q,  Ep)\ 


G(Q,EP) 


kT 


(2.2) 


The  term  [-U(Q,Ep)]  represents  the  potential  energy  of  the  chromophore  dipole.  It  is 
composed  of  two  energy  contributions.  The  first,  Uo,  is  from  the  intermolecular 
interactions  in  the  local  environment  of  the  chromophore.  The  second,  Ui,  arises  from 
the  chromophore/electric  poling  field  interactions.  Using  bold  type  to  represent  vector 
and  tensor  quantities,  Ui  is  expressed  as: 


U.-H.F-  I 


a*F*F 


(2.3) 


where  F  is  the  local  electric  field  at  the  chromophore.  A  typical  approximation  for 
(cos30)  involves  neglecting  the  chromophore-chromophore  intermolecular  electrostatic 
interactions  (Uo)  along  with  the  polarizability  term.  This  gives: 


Ui  =  |ifpEpCos0 


(2.4) 


where  p,  is  the  chromophore  dipole  moment,  Ep  is  the  magnitude  of  the  applied  electric 
poling  field,  and  fpEp  is  the  effective  poling  field  felt  by  the  chromophore.  The  term  fpEp 
is  equal  to  V/h,  where  V  is  the  voltage  of  the  applied  poling  field  and  h  is  the  film 
thickness.  Now  the  order  parameter  can  be  computed  according  to  statistical  mechanics: 


(cosx0)  =  Lx 


^fPEP 


.  ( —  uf  E  cosO^ 

Jexp  - — - -  cosx  0sin0d0 


r  r-iifnE.oosO^ 
fexp  - — - -  sin0d0 


where  the  integrals  are  evaluated  from  0  =  0  to  180°.  The  term  Lx(pfpEp/kT)  represents 
Langevin  functions,  which  are  related  to  the  averaged  values  of  the  Legendre 


polynomials  <Px(cos0)>.  Let  x  =  (pfpEp/kT)  to  get: 

(cos0>  =  Li(x)  =  coth(x)  -  (1/x)  =  <Pi(cos0))  (2.6) 

<cos20>  =  L2(x)  =  1  +  (2/x2)  -  (2/x)coth(x)  =  (1/3)[2(P2(cos0)>  +1]  (2.7) 

<cos30)  =  L3(x)  =  (1  +  6/x2)coth(x)  -  (3/x)(l  +  2/x2) 

=  (1/5)[2(P3(cos0)>  +  3<Pi(cos0>]  (2.8) 

The  Langevin  functions  can  be  expanded  in  a  power  series: 

Li(x)  =  (x/3)  -  (x3/45)  +  (2x5/945)  +  -  (2.9) 

L3(x)  =  Li(x)(l-6/x2)  -  (2/x)  -  (x/5)  -  (x3/105)  +■••  (2.10) 


Assume  pfpEp  «  kT  and  use  only  the  first  terms  of  the  power  series  to  get: 

(cos30)  =  x/5  =  pfpEp/5kT 


(2.11) 
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The  number  density,  N,  has  been  defined  by  workers  at  IBM,  Lockheed-Martin9, 
Akzo10  and  elsewhere  in  terms  of  the  weight  fraction  of  the  chromophore  in  the  polymer, 
w.  This  leads  to  the  relationship  between  r 33  and  p: 


2wpNAf,f.!E,  (  ^ 

15kTr)4  [mw> 


(2.12) 


where  NA  is  Avagadro’s  number,  p  is  the  density  of  the  nonlinear  polymer  system,  Mw  is 
the  molecular  weight  of  the  chromophore  and  T  is  the  poling  temperature.  This  suggests 
that  the  electro-optic  coefficient  should  increase  linearly  with  the  chromophore  number 
density  (loading).  However,  this  has  not  been  observed  with  experimentation.  For 
chromophores  with  large  polarizabilities  and  dipole  moments,  the  electro-optic 
coefficient  does  not  increase  linearly.  In  fact,  many  of  these  chromophores  exhibit 
maximum  peaks.  The  chromophore-chromophore  electrostatic  interactions  are 
responsible  for  the  deviation.  Figure  2.1  represents  this. 

Dalton  et.  al.  have  developed  a  theoretical  model  that  successfully  predicts  the 
trends  in  the  electro-optic  coefficient  as  the  loading  density  and  molecular  structure 
change.5  Dalton's  model  uses  ordering  and  disordering  forces  as  components  of  the 
acentric  ordering  parameter  (cos30).  With  so  many  intermolecular  electrostatic 
interactions  between  chromophores,  the  problem  is  a  many-body  problem.  Statistical 
mechanics  is  the  most  logical  way  to  approach  the  problem,  although  it  would  require 
solving  a  tremendous  number  of  rotational  matrices  that  relate  the  chromophores  to  each 
other  and  to  the  electric  poling  field.  The  integrals  of  equation  (2.5)  would  have  to  be 
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Fig.  2.1  Plot  of  theoretically  predicted  electro-optic 
coefficients  versus  number  density  for  the  FTC 
chromophore.  For  theoretical  calculations,  the  molecules 
are  treated  as  hard  spheres  or  prolate  ellipsoids.11 


solved  over  many  orientational  variables.  Instead,  a  coordinate  system  such  as  the  one 
shown  in  Figure  2.2  is  used  to  approximate  the  analytical  solution. 
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Ep  R(Q)  or  0 


Fig.  2.2  The  laboratory  reference  system  used  to  relate  chromophores  to  each  other  and  to  the 
electric  field.  R(Q)  is  the  rotation  matrix  relating  the  reference  chromophore  to  the  laboratory 
axis.  R(£20  is  the  rotation  matrix  for  the  interaction  between  the  reference  chromophore's 
principal  axis  and  the  director  axis  of  other  chromophores  that  contribute  to  the  field  felt  by  the 
reference  chromophore.11 

A  reference  chromophore  is  chosen  and  all  other  chromophores  are  arranged  in  a 
lattice  around  it,  where  the  chromophore  number  density  fixes  all  the  distances.  Now  the 
reference  chromophore  in  the  center  can  be  related  to  the  electric  poling  field  and  the 
other  chromophores.  By  choosing  a  director  axis,  all  the  chromophores  that  contribute  to 
the  field  felt  by  the  reference  chromophore  that  can  be  included  as  a  function  of  Qj.  An 
average  is  calculated  over  the  orientational  variables  related  to  the  other  chromophores 
per  Piekara,  which  is  analogous  to  London's  treatment.5,12  The  result  gives  an 
electrostatic  interaction  energy  between  chromophores,  U2,  depending  on  one  angle,  Q. 
The  equation  is: 


U2  =  -Wcos(Qi) 


(2.13) 


The  potential  function  used  in  the  ordering  parameter  (cos30)  is  then  given  by: 

U  =  Ui  +  U2  =  -pfpEpCosO  -  Wcos(0,)  (2.14) 

By  using  the  trigonometric  identity: 

cos0  =  cos0'cos0i  +  sin0'sin0icos((|>i  -  <)>')  (2.15) 

the  orientational  integrals  are  easier  to  solve.  Next,  the  chromophores  are  treated  as  hard 
spheres  to  give  an  expression  for  (cos  0)  as: 


(cos30) 


(2.16) 


where  L  is  the  Langevin  function  expressed  in  terms  of  W/kT  and  W  is  approximately 
the  intermolecular  electrostatic  interaction  energy  given  by  London.  As  the 
concentration  increases,  the  acentric  order  parameter  decreases.  Equation  (2.16)  verifies 
the  maxima  that  appear  in  the  plots  of  electro-optic  coefficient  versus  chromophore 
number  density.  From  Fig.  2.1  it  would  make  sense  to  produce  chromophores  with 
spherical  architecture,  which  has  been  done  successfully  by  Dalton  et.  al.5  proving  the 


theoretical  calculations  correct. 
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Chapter  3 

Design  and  Synthesis  of  Chromophores  Based  on  Reducing  Intermolecular 
Electrostatic  Interaction 

3.1  Introduction 

3.1.1  Motivation  for  Research 

One  of  the  biggest  challenges  in  nonlinear  optic  material  research  is  translating 
microscopic  to  macroscopic  optical  nonlinearity.  A  chromophore  can  be  thermally  and 
chemically  stable,  show  low  optical  loss,  display  high  pP  values  and  still  not  reflect  high 
r33  numbers  in  the  chromophore/polymer  composite  material.  As  mentioned  earlier, 
Dalton  and  coworkers  discovered  that  dipole-dipole  interactions  among  chromophores 
limit  the  degree  of  noncentrosymmetric  order  in  the  absence  of  chromophore 
modifications  that  sterically  inhibit  such  interactions.1,2  Last  year  Liakatas  et.  al?  and  Shi 
et.  al.4  reported  enhanced  V*  and  pP  values  based  on  bulky  donor  groups  and  modified 
bridges  to  reduce  intermolecular  electrostatic  interactions. 

3.2  Design  and  Synthesis  of  the  Chromophores 

3.2.1  Design  and  Synthesis  of  the  Triphenylamine  Donor 

There  are  three  reasons  for  choosing  this  particular  donor  structure.  The  first  is 
that  there  is  a  need  for  bulky  donor  groups  to  sterically  inhibit  interactions  with  other 
chromophores.  This  donor  should  also  not  have  any  intramolecular  H-bonding  that 
would  reduce  the  strength  of  the  donor.  The  second  reason  is  a  result  of  a  study  by 
Marder  and  Perry5  that  showed  an  increase  in  po  values,  although  slightly  red-shifted,  for 
di-aryl  amine  donors  over  di-alkyl  amine  donors,  both  attached  to  a  phenyl  group.  The 
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third  reason  is  a  result  of  research  done  at  EBM.6'10  Workers  there  found  an  increase  in 
thermal  stability  of  amine  donor  groups  by  switching  from  alkyl  amines  to  aryl  amines. 
The  donor  synthesis  was  started  with  a  variation  of  a  one-pot  palladium  catalyzed 
amination  published  by  Thayamanavan  et.  al.n  The  structure  and  synthetic  pathway  of 
the  donor  are  shown  in  Figures  3.1  and  3.2. 


Fig.  3.1  Structure  of  Triphenylamine  donor  in 
Wittip  salt  form 
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CH2PPh3Br 


Fig.  3.2  Synthetic  scheme  for  Triphenylamine  Donor 
in  Wittig  salt  form. 

3.2.2  Design  and  Synthesis  of  Chromophores  Based  on  Triphenylamine  Donor 

Each  chromophore  acceptor  utilized  the  2-dicyanomethylen-3-cyano-4,5,5- 
trimethyl-2,5dihydrofuran  acceptor  designed  at  Lockheed-Martin.  Marder12  16  and  Jen17 
20  have  also  shown  that  polyene  chain  and  thiophene-containing  bridges  give  higher  P 
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values.  Two  chromophores  were  produced  with  the  triphenylamine  donor  for 
comparison  to  two  other  similar  chromophores  containing  different  donors.  Their 
structures  are  shown  in  Figures  3.2  and  3.3.  The  synthesis  of  the  intermediates  followed 
literature  procedures  with  slight  modifications  as  shown  in  Figures  3.4, 3.5  and  3.6. 


MeO 


MeO 


Fig.  3.2  Structures  of  the  two  chromophores  with 
triphenylamine  donors  produced  for  this  study. 


46 


LDA,  THF,  0  °C 
N-formyl-piperidine 


HSCH2C02Et, 

- 

k2co3,  dmf 


LiOH,  THF 


Fig.  3.4  Synthetic  pathway  for  obtaining  the  dialdehyde 
version  of  thieno[3,2-b]thiophene  for  the  bridge. 


Fig.  3.5  Synthetic  pathway  for  obtaining  the  bridge  2,3- 
thiopheno[3,4-e]  1 ,4-dioxane-5,7-dicarbaldehyde 
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Fig.  3.6  Synthetic  pathway  for  obtaining  the  2- 
dicyanomethylen-3  -c  yano-4, 5 , 5 -trimethyl-2, 5  dihydrofuran 
acceptor. 

The  synthetic  schemes  for  the  chromophores  are  shown  in  Figures  3.7  and  3.8. 


Fig.  3.7  Synthetic  pathway  for  chromophore  5. 


48 


Fig.  3.8  Synthetic  pathway  for  chromophore  6. 

3.3  Experimental  Section 

3.3.1  General  Materials  and  Methods 

The  solvents  and  chemicals  used  were  purchased  from  Aldrich,  Acros  and  other 
chemical  companies.  They  were  analytical  reagent  grade  unless  otherwise  noted. 

Column  chromatography  was  performed  on  70-230  mesh  Lancaster  silica  gel  60. 
Conventional  ’H  NMR  spectra  were  recorded  using  a  Bruker  200  FT  NMR  spectrometer. 
The  chemical  shifts  are  referenced  to  tetramethylsilane  (TMS)  internal  standard. 
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3.3.2  Synthesis  of  Bis-(4-methoxyphenyl){4-[(triphenyIphosphonium  bromide) 
methyl]phenyl)amine  (1) 

Bis(4-methoxyphenyl)phenylamine  (2):  The  synthesis  of  the  triphenylamine 
donor  Wittig  salt  follows  the  schematic  in  Fig.  3.2.  To  a  solution  of  4.975g  (5.43  mmol) 
of  tris(dibenzylideneacetone)-dipalladium(0)  and  4.5 19g  (8.15mmol)  of  1,1’- 
bis(diphenylphosphine)-ferrocene  in  680mL  toluene  under  nitrogen  was  added  68.46mL 
(0.245mol)  of  4-bromoanisole  and  was  allowed  to  stir  for  25  minutes.  Then,  sodium 
tert-butoxide  (59.36g,  0.61 8mol)  and  aniline  (22.5mL,  0.236mol)  were  added  to  the 
solution  and  stirred  at  90°C  for  approximately  2  weeks.  Thin  layer  chromatography  was 
used  to  monitor  the  complete  formation  of  the  decoupled  product.  The  reaction  solution 
was  then  worked  up  with  brine  washings  (3x),  extracted  with  ether,  and  dried  over 
MgS04.  A  flash  column  of  5%  ethyl  acetate/95%  hexanes  gave  a  light  brown  solid.  The 
still  crude  product  was  purified  on  a  column  with  1%  ethyl  acetate/99%  hexanes  mobile 
phase  to  give  21.19g  of  white  solid.  lB  NMR  (CDC13):  83.87  (s,  6H),  86.87-7.08  (m, 
11H),  87.26  (d,  2H). 

4-[bis(4-methoxyphenyl)amino]benaldehyde  (3):  In  an  addition  funnel, 
0.34mL  (3.64mmol)  of  POCI3  was  added  dropwise  to  a  stirred  cooled  solution  at  0°C 
containing  0.76mL  (9.84mmol)  of  DMF  in  a  three-neck  flask  and  allowed  to  stir  for  one 
hour.  The  mixture  was  then  allowed  to  warm  to  room  temperature.  A  solution  of  2 
(l.Og,  3.28mmol)  in  1,2  dichloroethane  was  then  added  dropwise.  After  complete 
addition,  the  additional  funnel  was  replaced  with  a  condenser  and  the  solution  was 
heated  to  90-95°C  for  ~3  hours.  After  slight  cooling,  the  solution  was  added  dropwise  to 
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a  solution  of  NaHCCb.  The  crude  product  was  extracted  with  CH2CI2,  washed  3x  with 
NaHCC>3 ,  and  dried  over  Na2S04.  The  crude  product  was  purified  with  column 
chromatography  with  20%  ethyl  acetate/80%hexanes  as  the  mobile  phase  to  reveal  1 .04g 
of  viscous  bright  yellow  oil.  ]H  NMR  (CDCI3):  83.99  (s,  6H),  86.90  (d,  4H),  87.20  (d, 
4H),  87.61  (d,  2H),  89.79  (s,  1H). 

{4-[bis(4-methoxyphenyl)amino]phenyl}methan-l-ol  (4):  To  a  solution  of 
methanol,  14.9g  (45  mmol)  of  3  was  added  and  stirred.  To  a  prepared  solution  of  0.75g 
of  NaOH  in  2.5mL  H20  was  added  NaBH4  (0.85g,  22.5  mmol)  and  25mL  of  methanol. 
The  prepared  solution  was  added  to  the  stirred  solution  of  3  at  0°C  via  an  addition 
funnel.  The  solution  was  allowed  to  stir  at  room  temperature  overnight.  The  solution 
was  then  worked  up  with  brine  washings  (3x),  extracted  with  ether,  and  dried  over 
MgS04.  Removal  of  solvent  revealed  14.84g  of  dark  orange  oil.  The  product  was  used 
without  further  purification.  *H  NMR  (CDCI3):  83.95  (s,  6H),  84.73(s,  2H),  86.94  (d, 
4H)  87.07  (d,  4H),  87.20  (d,  2H),  87.31  (d,  2H). 

Bis-(4-methoxyphenyl){4-[(triphenylphosphonium 
bromide)methyl]phenyl)amine  (1):  A  solution  of  4  (14.84g,  0.044  mmol)  and 
triphenylphosphonium  hydrobromide  (13.7g,  0.040  mmol)  in  lOOmL  of  chloroform  was 
placed  on  an  azeotrope  distillation  apparatus  and  refluxed  to  remove  water  for  2-3  hours. 
Once  cooled  to  room  temperature,  the  chloroform  solution  was  concentrated  via  rotary 
evaporation.  The  product  was  precipitated  using  ether  and  27. 8g  were  isolated  by 
filtration.  The  'H  NMR  is  shown  in  Figure  3.9.  *H  NMR  (CDCI3):  83.80  (s,  6H),  85.22 
(d,  2H),  86.62  (d,  4H),  6.79  (d,  4H),  86.98  (d,  2H),  87.60-7.81  (m,  15H). 
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Fig.  3.9  ’H  NMR  spectrum  of  triphenylamine  donor  Wittig 
salt 


3.3.3  Synthesis  of  Chromophore  5 

Compound  10:  In  a  1L  round  bottom  flask,  6.67g  (10.09  mmol)  of  1  and  2.26g 
(20.14  mmol)  of  potassium  ferf-butoxide  were  added  to  about  300mL  of  methylene 
chloride.  The  solution  was  stirred  vigorously  for  about  two  minutes  at  room 
temperature.  Next,  83.70mg  (0.32  mmol)  of  18-crown-6  were  added  to  2g  of  8  and 
stirred  at  room  temperature  for  2-3  hours.  The  solution  was  filtered  over  Celite  and  roto- 
evaperated.  Column  chromatography  using  5%  ethyl  acetate  in  hexane  gave  l.lg  of 
product.  *H  NMR  (CDCI3):  83.91  (s,  6H),  54.42-4.59  (m,  2H),  86.91  (d,  4H),  57.09  (d, 
4H),  86.99  (d,  2H),  87.37  (d,  2H),  87.41  (d,  1H),  87.12  (d,  1H),  510.18  (s,  1H),  84.64  (q, 


2H). 
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Chromophore  5:  To  compound  10  just  enough  chloroform  to  dissolve  the 
product  was  added  to  a  small  vial  with  nearly  lg  of  10.  The  mixture  was  stirred  and 
allowed  to  reflux.  Then  about  2-3  drops  of  triethylamine  were  added,  and  the  solution 
continued  to  stir  for  about  30-45  minutes.  It  was  quenched  with  ammonium  chloride 
(2x)  and  washed  with  chloroform.  It  was  then  recrystallized  in  methanol  and  purified  by 
column  chromatography  using  3%  ethyl  acetate  in  hexanes  to  afford  0.6g  of  product.  !H 
NMR  (CDC13):  51.51  (s,  3H),  51.62  (s,  3H),  53.83  (s,  6H),  54.42-4.53  (m,  4H),  66.85  (d, 
1H),  66.90  (d,  4H),  57.46  (d,  1H),  57.13  (d,  1H). 

3.3.4  Synthesis  of  Chromophore  6 

Compound  11:  In  a  1L  round  bottom  flask,  5.47g  (8.28  mmol)  of  1  and  1.86g 
(16.56  mmol)  of  potassium  feri-butoxide  were  added  to  about  300mL  of  methylene 
chloride.  The  solution  was  stirred  vigorously  for  about  five  minutes  at  room 
temperature.  Next,  82.89mg  (0.31  mmol)  of  18-crown-6  were  added  to  2.5g  of  7  and 
stirred  at  room  temperature  for  2-3  hours.  The  solution  was  filtered  over  Celite  and  roto- 
evaperated.  Column  chromatography  using  methylene  chloride  gave  0.9g  of  product,  'll 
NMR  (CDCI3):  63.84  (s,  6H),  66.86  (d,  4H),  66.83  (d,  4H),  57.02  (d,  2H),  57.31  (d,  1H), 
57.46  (d,  2H),  58.02  (d,  1H),  59.98  (s,  1H). 

Chromophore  6:  To  compound  11  just  enough  chloroform  to  dissolve  the 
product  was  added  to  a  small  vial  with  0.86g  of  11.  The  mixture  was  stirred  and  allowed 
to  reflux.  Then  about  2-3  drops  of  triethylamine  were  added,  and  the  solution  continued 
to  stir  for  about  30-45  minutes.  It  was  quenched  with  ammonium  chloride  (2x)  and 
washed  with  chloroform.  It  was  then  recrystallized  in  methanol  and  purified  by  column 
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chromatography  using  methylene  chloride  to  afford  0.55g  of  product.  *H  NMR  (CDCI3): 
81.56  (s,  3H),  81.69  (s,  3H),  83.92  (s,  6H),  86.96  (d,  4H),  86.85  (d,  4H),  86.90  (d,  2H), 
88.24  (d,  1H),  87.93  (d,  1H). 

3.4  Conclusion 

The  bottom  line  for  these  syntheses  is  that  the  bulky  donor  unit  is  very  easy  to 
make  and  functionalize,  if  needed.  We  believe  that  it  will  offer  more  stability  to  the 
chromophore  in  both  physical  and  thermal  realms.  With  its  shape,  it  will  also  somewhat 
inhibit  intermolecular  electrostatic  interactions,  and  a  little  savings  between  the  energy  of 
the  chromophores  can  go  a  long  way  concerning  increased  electro-optic  coefficients. 
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